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BACKGROUND OF THE INVENTION 
The Field of the Invention 

[0012] This invention relates to the field of printed feature manufacturing, 
such as integrated circuit manufacturing. In particular, this invention relates to 
computing corrections using a proximity effects model at evaluation points to 
achieve improved agreement between a design layout and an actual printed feature. 

Description of Related Art 

[0013] To fabricate an integrated circuit (IC), engineers first use a logical 
electronic design automation (EDA) tool, also called a functional EDA tool, to 
create a schematic design, such as a schematic circuit design consisting of symbols 
representing individual devices coupled together to perform a certain function or set 
of functions. Such tools are available from CADENCE DESIGN SYSTEMS and 
from SYNOPSYS. The schematic design must be translated into a representation of 
the actual physical arrangement of materials upon completion, called a design 
layout. The design layout uses a physical EDA tool, such as those available from 
CADENCE and AVANT!. If materials must be arranged in multiple layers, as is 
typical for an IC, the design layout includes several design layers. 

[0014] After the arrangement of materials by layer is designed, a fabrication 
process is used to actually form material on each layer. That process includes a 
photo-lithographic process using a mask having opaque and transparent regions that 
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causes light to fall on photosensitive material in a desired pattern. After light is 
shined through the mask onto the photosensitive material, the light-sensitive material 
is subjected to a developing process to remove those portions exposed to light (or, 
alternatively, remove those portions not exposed to light). Etching, deposition, 
diffusion, or some other material altering process is then performed on the patterned 
layer until a particular material is formed with the desired pattern in the particular 
layer. The result of the process is some arrangement of material in each of one or 
more layers, here called printed features layers. 

[0015] Because of the characteristics of light in photolithographic 
equipment, and because of the properties of the material altering processes 
employed, the pattern of transparent and opaque areas on the mask is not the same as 
the pattern of materials on the printed layer. A mask design process is used, 
therefore, after the physical EDA process and before the fabrication process, to 
generate one or more mask layouts that differ from the design layers. When formed 
into one or more masks and used in a set of photolithographic processes and material 
altering processes, these mask layouts produce a printed features layer as close as 
possible to the design layer. 

[0016] The particular size of a feature that a design calls for is the feature's 
critical dimension. The resolution for the fabrication process corresponds to the 
minimum sized feature that the photolithographic process and the material processes 
can repeatably form on a substrate, such as a silicon wafer. As the critical 
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dimensions of the features on the design layers become smaller and approach the 
resolution of the fabrication process, the consistency between the mask and the 
printed features layer is significantly reduced. Specifically, it is observed that 
differences in the pattern of printed features from the mask depend upon the size and 
shape of the features on the mask and the proximity of the features to one another on 
the mask. Such differences are called proximity effects. 

[0017] Some causes of proximity effects are optical proximity effects, such 
as diffraction of light through the apertures of the optical systems and the patterns of 
circuits that resemble optical gratings. Optical proximity effects also include 
underexposure of concave corners (inside corners with interior angles greater than 
180 degrees) and overexposure of convex corners (outside corners with interior 
i t angles less than 1 80 degrees), where the polygon represents opaque: regions, and 

different exposures of small features compared to large features projected from the 
same mask. Other causes of proximity effects are non-optical proximity effects, 
such as sensitivity of feature size and shape to angle of attack from etching plasmas 
or deposition by sputtering during the material altering processes, which cause 
features to have shapes and sizes that have decayed from or accumulated onto their 
designed shapes and sizes. 

[0018] In attempts to compensate for proximity effects, the mask layouts can 
be modified. A proximity effects model is often developed and used to predict the 
location of edges in the printed layer, given a mask layout. 
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[0019] A proximity effects model is typically built for a particular suite of 
equipment and equipment settings assembled to perform the fabrication process. 
The model is often built by performing the fabrication process one or a few times 
with test patterns on one or more mask layouts, observing the actual features printed 
(for example with a scanning electron micrograph), and fitting a set of equations or 
matrix transformations that most nearly reproduce the locations of edges of features 
actually printed as output when the test pattern is provided as input. The output of 
the proximity effects model is often expressed as an optical intensity. Other 
proximity effects models provide a calibrated output that includes a variable 
intensity threshold as well as optical intensity. Some proximity effects models 
provide a calibrated output that indicates a printed edge location relative to an 
intensity threshold value as a spatial deviation. Thus, some model outputs can take 
on one or more values that vary practically continuously between some minimum 
value and some maximum value. 

[0020] Once a proximity effects model is produced, it can be used in the 
fabrication design process. For example, the proximity effects model is run with a 
proposed mask layout to produce a predicted printed features layer. The predicted 
printed features layer is compared to the design layer and the differences are 
analyzed. Based on the differences, modifications to the proposed mask layout are 
made. 

[0021] In conventional systems, the distances between the predicted edges 
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and the edges of the items on the design layout are used to correct the mask. For 
example, as shown in FIG. 1, the proximity model predicts that an original edge in 
the mask 1 10 aligned with a corresponding edge in the design layer, will result in a 
printed edge 120. The predicted printed edge 120 is displaced X units inward from 
the edge 1 10 as indicated by vector 132. For simplicity of illustration, the positions 
of the predicted printed features are shown on the same scale as the items in the 
design and fabrication layouts, both original and adjusted. The distance X units is 
measured on the common scale. The conventional system in this example modifies 
the mask layout by moving a segment corresponding to this edge 1 10 in the original 
mask layout outward, in the opposite direction, by a correction distance C equal to X 
distance units, as indicated by vector 134. The new edge in the adjusted mask is 
u ! denoted by segment 114. In other conventional systems, not shown, the correction 
vector C has a magnitude given by a constant proportion of X. 
[0022] FIG. 1 also shows that the difference between the predicted edge 1 20 
and the original fabrication layout edge 1 10 is computed from an evaluation point 
105 on a segment 112 of the edge 1 10. The segment 1 12 in the original fabrication 
layout is bounded by dissection points 102 and 104. The correction computed for 
the evaluation point 105, indicated by vector 134, is then applied to the whole 
segment 112, resulting in a new edge 1 14 for the adjusted mask layout. At least 
because of this change, the adjusted mask layout with edge 1 14 differs from the 
original fabrication layout with edge 110. The vertices 142 and 144 of the new edge 
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1 14 in the adjusted fabrication layout still correspond to the dissection points 102 
and 104 on the original fabrication layout edge. After making these corrections to all 
segments in the proposed, original mask layout, a final, adjusted mask layout is 
generated that is used to form a mask used in the fabrication process. 
[0023] A problem with the conventional technique is that changes in a 
printed features layer, such as a layer on a silicon wafer, may not be equivalent or 
linearly related to changes in position of an edge in a fabrication layout, such as a 
mask layout. As a consequence, a change in mask edge by X units, or proportional 
to X units, in general does not produce a movement of the printed feature edge 120 
by X units, even when mapped at the same scale. The relationship between changes 
in the mask and changes in the printed features is often referred to as the Mask Error 
Factor (MEF) or the Mask Error Enhancement Function (MEEF). For integrated 
circuited, the MEEF is often computed as the partial derivative of the wafer 
dimension with respect to the mask dimension on the common scale. 
[0024] The deviation of the MEEF from a constant value often prevents a 
mask corrected in the conventional way from fabricating printed features that meet 
the spatial tolerance required by a manufacturing specification, i.e., the spec 
tolerance. 

[0025] In this field, therefore, new ways are needed to compute corrections 
at evaluation points, that account for non-equivalence and non-linearity in the 
MEEF, so that the printed edges are guaranteed to lie within a spec tolerance of the 
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design layer whenever such is possible for a given suite of fabrication equipment. 
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SUMMARY OF THE INVENTION 

[0026] According to techniques of the present invention, an optimal bias is 
computed for a segment in a mask layout using a proximity effects model to account 
for non-equivalence and non-linearity between changes in mask layouts and changes 
in printed features. As used here, a bias is a measure of the distance between a 
segment of an edge in an adjusted fabrication layout, like a mask layout, from a 
corresponding edge in the original fabrication layout. The edge in the original 
fabrication layout corresponds to an edge in a design layout. The optimal bias is 
selected to give an adjusted mask layout that produces an edge in the printed 
features layer in closer agreement with a corresponding edge in the design layer than 
is available using conventional procedures. 

[0027] * In particular, in several embodiments, the proximity effects model is 
run for more than one value of bias, and the optimal bias is deduced from the output 
of these two or more model runs. Several embodiments of these techniques include 
elements that speed the computation of the proximity effects model and that select 
the two or more biases so that the computations can be applied over an entire layout 
iteratively, and still rapidly converge on optimal biases for the segments in the 
layout. 

[0028] In other embodiments, an inverse proximity effects model is run. The 
positions of edges in segments in the design are input to the inverse model and 
amplitudes related to an edge in the fabrication layout are output at several locations. 
The optimal bias is deduced from the output of the inverse model. 
[0029] In one aspect of these techniques, a technique of correcting for 
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proximity effects a fabrication layout corresponding to an original fabrication layout 
for a printed features layer includes executing a proximity effects model to produce 
an initial output. The initial output is based on a first position for a segment in a 
fabrication. layout. The first position is displaced from a corresponding original edge 
in the original fabrication layout by a distance equal to an initial bias. The model is 
also executed to produce a second output based on a second position for the segment. 
The second position is displaced from the corresponding original edge by a distance 
equal to a second bias. An optimal bias for the segment is determined based on the 
initial output and the second output. The segment is displaced in the fabrication 
layout from the corresponding edge based on the optimal bias. 
[0030] In some embodiments of this aspect, the method further comprises 
determining a first correction based on the initial output associated with the initial 
< bias. A second correction is determined based on the second output and associated 
with the second bias. Determining the optimal bias includes interpolating to a third 
bias associated with a target correction, and setting the optimal bias based on the 
third bias. Each correction is related to a difference between a first value of a 
property at a predicted printed edge and a second value of the property on the 
original edge. 

[0031] In some embodiments of this aspect, executing the model and 
determining the optimal bias on the segment are repeated for a number of iterations 
until a termination criterion is satisfied. 

[0032] In some of these embodiments, executing the model and determining 
the optimal bias on the segment are repeated for a plurality of segments in the 
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fabrication layout during each iteration of the number of iterations. 
[0033] In another aspect of these techniques, a technique of correcting for 
proximity effects a fabrication layout corresponding to an original fabrication layout 
for a printed features layer includes executing an inverse proximity effects model for 
a segment of an edge in the original fabrication layer. A correction is determined for 
the segment based on a difference between an output from the inverse proximity 
effects model and an initial position for the segment in the fabrication layout. 
[0034] In another aspect of these techniques, a technique of correcting for 
proximity effects a fabrication layout corresponding to an original fabrication layout 
for a printed features layer includes building an inverse proximity effects model. 
The inverse proximity effects model is built using information from a test mask with 
a plurality of polygons of uniform light transmission properties and a printed 
features layer fabricated using the test mask. The inverse proximity effects model 
accepts as input positions of edges in the printed features layer and produces an 
output related to a position of an edge on the test mask. The inverse proximity 
effects model is executed for a segment of an edge in the original fabrication layer. 
A correction is determined for the segment based on a difference between an output 
from the inverse proximity effects model and an initial position for the segment in 
the fabrication layout. 

[0035] In another aspect of these techniques, a technique for building an 
inverse proximity effects model includes using information from a test mask with a 
plurality of polygons of uniform light transmission properties and a printed features 
layer fabricated using the test mask. The inverse proximity effects model accepts as 
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input positions of edges in the printed features layer and produces an output related 
to a position of an edge on the test mask. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0036] The present invention is illustrated by way of example, and not by 
way of limitation, in the figures of the accompanying drawings and in which like 
reference numerals refer to similar elements and in which: 

[0037] FIG. 1 is a plan view of an edge on a mask and the resulting edge of a 
printed features layer predicted by a proximity effects model, showing that the 
difference is used to correct the mask according to conventional methods. 
[0038] FIG. 2 is a flow diagram showing the relationship among sub 
processes in a printed features layer design and fabrication process, according to 
embodiments of the present invention. 

[0039] FIG. 3 A is a plan view of a portion of a design layer showing 
orientation of a section for producing profiles of amplitude output by model runs. 
[0040] FIG. 3B is a graph plotting profiles along the section of FIG. 3 A for 
several different values of bias. 

[0041] FIG. 3C is a plan view of an adjustment to a fabrication layout. 
[0042] FIG. 4A is a flow diagram showing the iterative adjustment of a 
fabrication layout using a proximity effects model. 

[0043] FIG. 4B is a flow diagram showing details of step 402 in FIG. 4A, 
according to one embodiment. 

[0044] FIG. 4C is a flow diagram showing details of step 450 in FIG. 4 A, 
according to an embodiment. 
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[0045] FIG. 4D is a flow diagram showing details of step 408 in FIG. 4A, 
according to an embodiment. 

[0046] FIG. 4E is a flow diagram showing details of step 410 in FIG. 4A, 
according to an embodiment. 

[0047] FIG. 4F is a flow diagram showing additional details of step 410 in 
FIG. 4A, according to an embodiment. 

[0048] FIG. 4G is a flow diagram showing additional details of step 410 in 
FIG. 4A, according to an embodiment. 

[0049] FIG. 4H is a flow diagram showing details of step 430 in FIG. 4A, 
according to an embodiment. 

[0050] FIG. 5 A is a flow diagram showing details of step 412 in FIG. 4B, 
according to an embodiment. 

[0051] FIG. 5B is a flow diagram showing details of step 41 1 in FIG. 4E, 
according to an embodiment. 

[0052] FIG. 6 is a flow diagram for building an inverse proximity effects 
model, according to an embodiment of another aspect. 
[0053] FIG. 7 is a flow diagram showing the iterative adjustment of a 
fabrication layout using an inverse proximity effects model, according to one 
embodiment of another aspect. 

[0054] FIG. 8 is a block model of a computer system configured according 
to one embodiment. 
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[0055] FIG. 9 is a plan view of an adjustment to corner segments in a 
fabrication layout. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0056] Techniques are described for correcting fabrication layouts for 
proximity effects for the fabrication of printed features layers, such as in integrated 
circuits. 

FUNCTIONAL OVERVIEW 

[0057] FIG. 2 shows the processes and layouts produced according to the 
present techniques for designing and fabricating printed features layers 249. The 
conventional processes are modified to include new techniques for modeling 
proximity effects and adjusting a proposed mask layout, as represented in FIG. 2 by 
process 263. 

[0058] The conventional processes include the Functional EDA process 210 
that produces the schematic diagram 215. The physical EDA process 220 converts 
the schematic diagram to a design layout made up of one or more design layers 225. 
After mask layouts 235 are produced, the conventional processes employ a 
fabrication process 240 to produce the printed features layer 249 using the mask 
layouts 235. The printed features layer 249 may be a printed circuit or the mask used 
to produce the printed circuit. In the former case, the fabrication process includes 
one process 243 for forming the mask and a second process 245 to produce the layer 
of the integrated circuit using the mask. If the printed features layer is the mask, step 
245 is skipped. The device 250, such as an integrated circuit, comprises one or more 
printed features layers. 



17 



NTI-019-5-1D 



PATENT 



[0059] According to techniques of the present invention, the fabrication 
design process 260 includes a new process 263 for modeling proximity effects and 
adjusting the fabrication layout. 

[0060] The modeling is done at evaluation points using as input the 
distribution of polygons in a fabrication layout, such as the original proposed mask 
layout 261 . Each evaluation point is placed on a segment of an edge on the original 
proposed layout 261 which corresponds to one of the design layers, e.g., 225a. A 
segment is defined by a pair of dissection points associated with the evaluation point 
on the edge in the original fabrication layout. The dissection points and evaluation 
points are selected in step 262 based on the polygons in the fabrication layout that 
are associated with the polygons in the design layout. 

[0061] In a preferred embodiment, the determination of dissection points and 
evaluation points on the fabrication layout proceed as described in the related 
applications, cited above. Each evaluation point and dissection point so determined 
is related to. corresponding points in the original fabrication layout. For example, in 
FIG. 1, the dissection points 102 and 104 on edge 1 10 of the original fabrication 
layout correspond to the vertices 142 and 144, respectively, of the segment 1 14 in 
the adjusted fabrication layout. 

[0062] The proximity effects are modeled in process 264. Based on an 
analysis of the model results and a comparison with the original fabrication layout 
corresponding to the design layer, as described in detail below, an optimal bias is 
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computed in this process. The optimal bias is then applied to an entire segment in 
the original fabrication layout corresponding to the segment in the design layer to 
form a corresponding segment in the adjusted fabrication layout. The optimal bias 
value expresses the distance from the position of the segment in the original 
fabrication layout to the position of the corrected segment in the adjusted fabrication 
layout with a positive value for an outward adjustment. As a result of the 
adjustments, the polygons on the adjusted mask layout may differ from the polygons 
in the original proposed mask layout. The final mask layouts 235 are then based on 
the adjusted mask layouts produced by the adjustment process 263. In an iterative 
implementation, the result of the adjustment process becomes a new (non-original) 
proposed layout that is again subjected to the adjustment process until some 
termination criterion is satisfied. 

ADJUSTING A SINGLE SEGMENT 

[0063] To explain several embodiments of the process 264, it is first 
described in one embodiment in which it is applied with a conventional proximity 
effects model to a single segment after the dissection points and the evaluation 
points have been selected for the layout. In the embodiment described here, 
considerations of execution time are not included, so this embodiment is not limited 
in the number of points where the proximity effects model may be run. This 
embodiment is described with respect to FIG.s 3 A, 3B and 3C. 
[0064] In the following, the term amplitude refers to a scalar output from a 
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proximity effects model. The term amplitude includes the optical intensity, the 
spatial deviation, the variable threshold, any other value output from the model, or 
any combination of these, such as a difference between an intensity and a variable 
intensity threshold, depending on the proximity effects model. In this terminology, 
an amplitude threshold is a value of the amplitude at a location where an edge is 
predicted to print. 

[0065] FIG. 3 A is a plan view of a portion of an original fabrication layout 
showing two opaque areas represented by polygons 302, 304. On the common scale 
used in this example, the arrows 307a, 307b and 307c represent a distance of 200 
nanometers (nm). Polygon 304 has an edge 303 with a segment defined by 
dissection points 305a and 305b. On this segment is evaluation point 306. Section 
308 intersects the segment at the evaluation point and is perpendicular to it. The 
amplitudes output by a proximity effects model are shown as profiles along section. 
308 from point 309a at 100 nm from the beginning of the section to point 309b at 
500 nm from the beginning. 

[0066] FIG. 3B is a graph plotting profiles along section 308 from FIG. 3 A 
for several different values of bias applied to the corresponding segment in the 
fabrication layout. The horizontal axis 314 gives the distance from the beginning of 
section 308. The vertical axis 312 gives the value of the amplitudes output by the 
proximity effects model. This graph plots 1 1 profiles of amplitudes output by a 
proximity effects model as curves 325 to 335. 
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[0067] Each profile represents the model output along section 308 for a 
different placement of the segment in an adjusted fabrication layout. The symbol in 
the curve is related to the placement of the segment according to the legend to the 
right of the graph. In the legend, the symbols are listed in column 316 and the 
associated placement of the segment in the fabrication layout is listed in column 
317. 

[0068] In this example, when the corresponding segment in the adjusted 
fabrication layout is placed at 400 nm, the segment coincides with the segment in the 
original fabrication layout, and the bias is zero. Thus curve 330 represents the 
profile output by the proximity effects model when the segment in the fabrication 
layout is at 400 nm; i.e., when the segment in the fabrication layout coincides with 
the original fabrication layout. In this description a negative bias refers to a 
displacement of a segment inside the corresponding edge in the original fabrication 
layout. In FIG. 3B a negative bias moves the segment to the right where the 
distance along section 308 increases. A positive bias refers to a displacement of a 
segment outside the corresponding edge, toward the left in FIG. 3B where the 
distance along section 308 decreases. 

[0069] To illustrate this method and distinguish conventional methods, the 
model is described as if it is a calibrated model with variable intensity threshold. 
For this illustration, the amplitude is taken to represent the optical intensity. For 
purposes of illustration, the value of the variable intensity threshold is taken to be 
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0.4 amplitude units near the evaluation point, and the amplitude threshold is set 
equal to the variable intensity threshold. Similar results are obtained if the 
amplitude is defined as the difference between the intensity and the variable 
intensity threshold, and a constant amplitude threshold is taken to be zero, and the 
curves in FIG. 3B are modified accordingly. 

[0070] Vertical line 315 marks the 400 nm distance where the segment lies 
in the original fabrication layout. Line 315 marks the nominal edge position, also 
called the target edge position for this segment. As can be seen from FIG. 3B, at a 
distance of 400 nm, curve 330 has a value substantially less than 0.4 amplitude 
units. Thus, the model predicts an edge is not printed at 400 nm in the printed 
features layer. Instead, the model predicts the edge is printed where the curve 330 
crosses the horizontal line at 0.4 amplitude units. The curve 330 crosses the 
horizontal line of 0.4 amplitude units at 378.9 nm, a distance 21.1 nm to the left of 
the desired target edge location in FIG. 3B. 

[0071] According to the conventional methods, a bias of -2 1 . 1 nm would be 
applied to move this segment in the fabrication layout to the right, i.e., to 421.1 nm. 
The hope is that moving the segment in the fabrication layout to the right by 21 .1 
nm will cause the edge in the printed layer to move to the right by the same amount. 
However, FIG. 3 B. shows that if the segment in the adjusted fabrication layout is 
moved to 420 nm, the profile will follow curve 332, which has a value at 400 nm 
that is greater than 0.4 amplitude units. Thus moving the segment in the fabrication 
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layout by even 20 nm is too far. The bias of -21.1 nm is too much because changes 
in the printed features layer are not equivalent to or linearly related to changes in the 
mask layout. 

[0072] According to this embodiment of the present invention, the proximity 
effects model is run for several different experimental values of bias for the segment 
in the adjusted fabrication layout. For example, the proximity effects model can be 
run for eleven different experimental values of the bias, from an experimental bias 
of 50 nm represented by curve 325 to an experimental bias of -50 nm represented by 
curve 335. An estimate of the optimal bias for this segment is obtained by 
interpolating between two of the profiles that straddle the intensity threshold value 
of 0.4 amplitude units at the target edge location, 400 nm. As can be seen in FIG. 
3B, curve 331 has a value at 400 nm that is less than the amplitude threshold value 
of 0.4 amplitude units; while the next profile, curve 332, has a value at 400 nm that 
is greater than the amplitude threshold value. In this example, an amplitude of 0.4 at 
400 nm is 65 percent of the way from curve 331 to curve 332 at 400 nm. According 
to this embodiment, a bias corresponding to a segment position 65 percent of the 
way from 410 nm to 420 nm is selected as the optimal bias. That is, the optimal bias 
is -16.5 nm according to this embodiment. 

[0073] Curve 340 in FIG. 3B plots the amplitude at a position of 400 nm, the 
nominal or target edge position, as a function of the position of the segment in the 
adjusted fabrication layout. Curve 340 intersects the amplitude threshold value of 
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0.4 amplitude units at 416.5 nm. The interpolating between profiles to the amplitude 
threshold value, described in the preceding paragraph, is equivalent to interpolating 
along curve 340 to the amplitude threshold value. In either case, applying a bias of 
[0074] -16.5 nm to the segment in the fabrication layout, rather than the - 
21.1 nm of the conventional methods, provides an edge in the printed features layer 
that is predicted to be closer to the target edge in the design layer, according to this 
proximity effects model. 

[0075] For a calibrated model that uses a threshold and a spatial deviation, a 
curve (not shown) that plots spatial deviation at the evaluation point for several 
different experimental biases can be used in lieu of curve 340. Then the 
interpolation is carried out to find a bias that yields a spatial deviation of zero. The 
. optimal bias is then based on this interpolated bias. 
[0076] FIG. 3C shows the relationship between the optimal bias computed 
according to the present invention and the bias applied according to the conventional 
methods. Vector 132 connects the evaluation point 105 to the predicted printed edge 
120 when no bias is applied. Vector 134 represents the bias according to the 
conventional methods. It is equal or proportional to vector 132 but in the opposite 
direction. For example, if vector 132 represents a distance of 21.1 nm to the left in 
FIG. 3B, then vector 1 34 represents a distance of 2 1 . 1 nm to the right, i.e., a bias of 
-21.1. In contrast, vector 394 represents the optimal bias applied to obtain adjusted 
segment 396 in the fabrication layout according to this embodiment of the present 
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invention. The magnitude of vector 394 is different than the magnitude of vector 
134. For example, as shown in FIG. 3C, vector 394 is smaller in magnitude than 
vector 132, with a magnitude of 16.5 nm rather than 21.1 nm. The optimal bias 
computed in this embodiment, and shown in FIG. 3C, accounts for the non- 
equivalence and non-linearity of the proximity effects captured by the proximity 
effects model. The conventional bias does not. 

REDUCING THE COMPUTATIONAL LOAD 

[0077] The techniques of the present invention include embodiments that use 
fewer computational resources than used to produce several full profiles, such as 
graphed in FIG. 3B. There, 1 1 profiles were generated, one for each trial value of 
bias from -50 to +50 in increments of 10 nm. For each profile the model was run 41 
times, one run for each location on section 308 from 100 to 500 nm in increments of 
10 nm. Thus 451 model runs were required to produce FIG. 3B. It is not necessary 
to run the model that many times to practice the techniques of the present invention. 
[0078] In one embodiment, each profile has a single value produced by 
running the proximity effects model at the evaluation point. Since the evaluation 
point is at the target position for the edge a curve for interpolation is generated 
directly. For example, curve 340 is produced simply by running the proximity 
effects model at the evaluation point for eleven experimental biases. Similarly, a 
plot of spatial deviation at the evaluation point can be produced by running that 
proximity effects model for the eleven experimental biases. The optimal bias is 
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obtained by interpolating along the resulting curve to the amplitude threshold value, 
either the intensity threshold value or zero spatial deviation, respectively. Thus the 
curve to be interpolated can be produced with only 1 1 model runs. This amounts to 
a savings of 440 model runs. 

[0079] In another embodiment, the model is run at the evaluation point only 
for two experimental values of the bias. For example, the two experimental biases 
are selected as zero and 20 nm, corresponding to a single value from each of the two 
profiles represented by curves 330 and 332. This embodiment then requires only 2 
model runs, a savings of 449 model runs. However, if the two experimental biases 
are too far apart, or they do not straddle the amplitude threshold value, such as zero 
spatial deviation, then the optimal bias obtained by interpolation equations may not 
be very precise. If this embodiment is used; it is important to select the two 
experimental biases carefully so that the optimal bias is precise within a 
specification tolerance whenever possible. Details on selecting the experimental 
bias values for a particular embodiment are described in more detail below. 
[0080] Other computational savings can be obtained if less then all the 
kernel functions in the proximity effects model can be used to produce the 
amplitudes, or if the proximity effects model can be run in consideration of only 
changes in the fabrication layout, or in consideration of only the nearest neighbors in 
the layout, or any combination of these savings. Techniques to realize these savings, 
as well as other advantages of the present invention, are illustrated with respect to a 
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particular embodiment of the invention described below. 

EXAMPLE EMBODIMENT WITH PROXIMITY EFFECTS MODEL 
[0081] FIG. 4A is a flow diagram showing an iterative adjustment of a 
fabrication layout using a proximity effects model, according to one embodiment of 
adjustment process 263 from FIG. 2. 

[0082] In step 402, parameters used by the techniques are initialized. More 
details on the parameters and their initialization are given below. In step 404, the 
current iteration begins. A segment of the fabrication layout is visited only once per 
iteration. More details on the steps included within step 404 are given below. In 
step 406, one of the segments not yet visited in this iteration is made the current 
segment. In step 408, it is determined whether the current segment should be 
skipped. If the current segment is skipped, control passes to step 438 to determine 
whether another segment remains to be visited in the layout. More details on the 
steps included within step 408 are given below. If the current segment is not 
skipped, control passes to step 410 in which the proximity effects model output is 
used for at least two different values of bias for the current segment. Typically, this 
involves executing once or twice a computer routine implementing the proximity 
effects model. Also in step 410, an optimal bias is determined based on the outputs 
from the two or more model outputs. Extensive details on the steps within step 410 
are given below. In step 430, the layout is adjusted if necessary by moving the 
current segment to a position separated from the corresponding segment in the 
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original fabrication layout by the optimal bias. More details on step 430 are given 
below. 

[0083] After it is determined in step 438 that another segment remains to be 
visited within the layout, control passes back to step 406 to make the next segment 
the current segment. If it is determined in step 438, that no other segment remains to 
be visited, then control passes to step 440 to determine whether all segments in the 
layout require corrections less than the specification tolerance. t If so, then no further 
changes to the biases are required, and the process ends in step 490. If at least one 
segment requires a correction greater than the specification tolerance, then control 
passes to step 450 to increment to the next iteration, and Visit at least one of the 
segments again. 

INITIALIZATION OF ADJUSTMENT PARAMETERS 

[0084] FIG. 4B is a flow diagram showing details for embodiment 402a of 
the initialize parameters step 402 in FIG. 4 A, 

[0085 ] In step 400, an initial value for the experimental additional bias is set. 
The experimental additional bias is a distance added to a first experimental bias to 
obtain a second experimental bias. The techniques of the present invention require 
at least two experimental biases for a segment in the fabrication layout, one of which 
could be zero, as explained above. As also explained above, the precision of the 
optimal bias obtained depends at least in part on the difference in the two 
experimental biases. The experimental additional bias parameter determines the 



28 



NTI-019-5-1D 



PATENT 



difference in the two experimental biases. This parameter is also called the 
experimental bias step size. 

[0086] In this embodiment, the experimental additional bias parameter is 
allowed to change with each iteration, so that step 400 sets a value for the first 
iteration, at least. In this embodiment, the experimental additional bias decreases 
from a maximum experimental additional bias on the first iteration toward a 
minimum experimental additional bias on subsequent iterations. In step 400, both 
the maximum experimental additional bias and the minimum experimental 
additional bias are set. 

[0087] The experimental additional bias is allowed to change because, on the 
first iteration, the amount of correction required is expected to be greatest. Also, a 
large value of the experimental additional bias is desirable to increase the likelihood 
that the two experimental biases will produce predicted edges that straddle the target 
edge position. Straddling the target edge insures that the bias associated with the 
threshold value is obtained by an interpolation rather than an extrapolation 
computation. Although the mathematical formula for extrapolation is often the same 
as for interpolation, in general it is not desirable to perform an extrapolation 
computation because it can lead to uncontrolled behavior. For example, the 
extrapolated bias is not bounded by the two experimental biases, and can become 
arbitrarily large. Because other iterations will follow, it is not necessary that the first 
few iterations produce the most accurate value for the bias. It is more important that 
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the first few iterations straddle the threshold value. To improve the precision during 
the later iterations, the experimental additional bias is decreased in the later 
iterations. The rate at which the experimental additional bias decreases with 
iteration is determined in this embodiment by a fading factor, fl, as is described in 
more detail below with respect to Equation 1 . In this embodiment, a value for the 
fading factor value fl is set during step 400. 

[0088] In step 40 1 , an initial value for the halo distance is set. The halo 
distance determines how far another feature in the fabrication layout may be from 
the evaluation point and still be used as input to the proximity effects model run at 
the evaluation point. In this embodiment, the halo distance is allowed to change 
based on iteration number. That is, the halo distance decreases from a maximum 
halo distance on the first iteration toward a minimum halo distance on subsequent 
iterations. In step 401, both the maximum halo distance and the minimum halo 
distance are set. The halo distance is allowed to decrease because the changes to the 
biases for the segments are expected to become smaller and fewer on later iterations. 
The smaller bias changes means that the segments are moving less on the later 
iterations. A smaller movement of a segment is likely to impact only other segments 
that are close by. Thus very little is added to the amplitude at a point by the more 
distant segments, when the segments have moved only a little. By decreasing the 
halo distance, fewer of the distant segments are involved in the computation. This 
decreases the computational effort to run the model without giving up significant 
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precision. The rate at which the halo distance decreases with iteration is determined 
in this embodiment by a fading factor, f2, as is described in more detail below with 
respect to Equation 2. In this embodiment, a value for the fading factor value f2 is 
set during step 401. 

[0089] In step 403, an initial value for the number of kernel functions used in 
running the proximity effects model is set. The number of kernel functions used 
determines the accuracy of the proximity effects model. The proximity effects model 
is built with a set of M kernel functions, where M is an integer. The number of 
kernel functions needed to achieve a given accuracy, such as the specification 
tolerance, is expected to decrease with each iteration because the geometries 
involved get smaller with each iteration. The geometries become smaller with each 
iteration because changes in segment biases get smaller with each iteration for the 
reasons given above. Also, the segments themselves may get smaller as new 
dissection points may be added to some segments on some iterations. As the 
geometries get smaller, the contributions of some kernel functions to the total model 
amplitude become negligible compared to the total. These kernel functions can be 
dropped. In this embodiment, a minimum number N of kernel functions is set for 
use after a large number of iterations, where N is an integer less than M. Decreasing 
the number of kernel functions used decreases the computational effort required to 
run the proximity effects model. The rate at which the number of kernel functions N 
decreases with iteration is determined in this embodiment by a fading factor, £3, as is 
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described in more detail below with respect to Equation 3. In this embodiment, a 
value for the fading factor value f3 is set during step 403. 
[0090] In step 405, a pre-bias is set for each segment. With a priori 
knowledge, a non-zero pre-bias can be used on at least some segments to get those 
segments closer to their ideal position before the first iteration. This reduces the 
number of iterations required to satisfy a termination criterion that depends on the 
changes in biases between successive iterations. 

[0091] In step 407, a value for the maximum bias is set. The maximum bias 
is a limit placed on an optimal bias, so that the optimal bias never becomes larger 
than the maximum bias. The maximum bias can be set as a fixed value, or it can be 
set as a fixed fraction of the distance separating neighboring edges in the fabrication 
layout. For example, referring to Figure 3A, the bias for segments on edge 303 of 
polygon 304 might reasonably be set to one third of the distance separating edge 303 
and the closest edge of polygon 302. In this case the fixed fraction for the maximum 
is set to substantially 1/3. In one embodiment, a maximum for positive bias is set 
separately from another maximum for a negative bias during step 407. 
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UPDATING THE ADJUSTMENT PARAMETERS EACH ITERATION 
[0092] FIG. 4C is a flow diagram showing details of one embodiment of the 
increment iteration step 450 of FIG. 4A. In this embodiment 450a, the experimental 
additional bias, the halo distance, and the number N of kernel functions are updated 
in each iteration. 

[0093] In step 45 1 , the experimental additional bias is updated. Any method 
that changes the experimental additional bias can be used in other embodiments. In 
the preferred embodiments, the experimental additional bias is decreased in later 
iterations. According to one embodiment, the experimental additional bias is 
decreased as follows: 

di = d m i n + (d max - d min ) / fl 1 (1) 
where I is the iteration number and begins with a value of zero, incrementing by one 
with each iteration; di is the experimental additional bias on iteration I; d m j n is the 
minimum experimental additional bias; d ma x is the maximum experimental 
additional bias; and fl is a first fading factor having a value no smaller than 1. A 
value for the first fading factor is chosen in order to control the speed with which the 
experimental additional bias decreases to the minimum experimental additional bias. 
The larger the first fading factor, the greater the speed with which the experimental 
additional bias decreases. In this embodiment, a value of fl is set in step 400 of 
FIG. 4B. 
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[0094] In step 453, the halo distance is updated. Any method that changes 
the halo distance can be used in other embodiments. In the preferred embodiments, 
the halo distance is decreased in later iterations. According to one embodiment, the 
halo distance is decreased as follows: 

hi = h m in + (h max - h min ) / ft 1 (2) 
where I is the iteration number; hi is the halo distance on iteration I; h m j n is the 
minimum halo distance; h max is the maximum halo distance; and f2 is a second 
fading factor having a value greater than 1 . A value for the second fading factor is 
chosen in order to control the speed with which the halo distance decreases to the 
minimum halo distance. The larger the second fading factor, the greater the speed 
with which the halo distance decreases. In. this embodiment, a value of £2 is set in 
step 401 of FIG. 4B. 

[0095] In step 455, the number N of kernel functions is updated. Any 
method that changes the number N of kernel functions can be used in other 
embodiments. In the preferred embodiments, the number N of kernel functions is 
decreases based on iteration umber I. According to one embodiment, the number N 
of kernel functions is decreased as follows: 

N, = INTEGER [ N min + (M - N min ) / fi 1 ] (3) 
where I is the iteration number; Nj is the number N of kernel functions on iteration I; 
N m j n is the minimum number N of kernel functions; M is the maximum number of 
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kernel functions; and D is a third fading factor having a value greater than 1 . The 
third fading factor is chosen in order to control the speed with which the number N 
of kernel functions decreases to the minimum number N of kernel functions. The 
larger the third fading factor, the greater the speed with which the number N of 
kernel functions decreases. In this embodiment, a value of f3 is set in step 403 of 
FIG. 4B. INTEGER[ ] indicates an operation which converts a value inside the 
brackets to an integer value, such as the least greater integer or the greatest lesser 
integer. This operation is useful because an integer number of kernel functions must 
be used. 

[0096] By the time control passes to step 450a, the first iteration has been 
completed and at least one segment has had the proximity effects model run in the 
vicinity of its evaluation point. If the model output includes a pre-image value, this 
pre-image associated with the evaluation point value is stored in step 457 in this 
embodiment. In many implementations of proximity effects models, an intensity is 
computed by first convolving the kernel functions with geometric characteristics of 
the edges in the fabrication layout relative to the evaluation point to produce a 
quantity called a pre-image. This portion of the process is linear and can be reused 
in different iterations and experimental biases. The pre-image values are used in a 
non-linear sum of squares to produce intensity. Since optical intensity on an two- 
dimensional array yields an optical black and white image, the value computed 
before intensity is a reasonably called a "pre-image" value. Characteristics of the 
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intensity or pre-image or both are then used in a calibrated model to predict where 
the edge will print, for example, to compute either a variable threshold or a spatial 
deviation. In proximity effects models that provide a pre-image value, the pre- 
image value can be used to speed up the computation at the evaluation point on the 
next iteration or with the next experimental bias, as described in more detail below. 
In any case, if available, the pre-image associated with the evaluation point is stored 
in step 457 according to this embodiment. In embodiments in which one or more 
supplemental evaluation points are associated with the evaluation point, the pre- 
image is stored for the one or more supplemental evaluation points. Supplemental 
evaluation points are described in more detail below. 

[0097] In another embodiment, new dissection points and evaluation points 
are determined on segments that- have been displaced during the current iteration. 
For example, in FIG. 2, process 262 for obtaining dissection and evaluation points is 
repeated on every iteration of the adjustment process 263. In this case, for example, 
a displaced segment in the adjusted fabrication layout, like segment 196 in FIG. 3C, 
may be close enough to a vertex of another edge to warrant further dissection. The 
new dissection points (not shown) determined as if the edge were at the position of 
segment 396 are then associated with new dissection points (not shown) on the 
original fabrication layout edge 1 10 to define new segments (not shown) on that 
original edge 110. 

SKIPPING SEGMENTS 
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[0098] FIG. 4D is a flow diagram showing details for embodiment 408a of 
the skip segment step 408 in FIG. 4A. In step 481, other segments within the 
current halo distance are reviewed. During this step, information is obtained that 
indicates whether a particular segment has had its bias changed during the last 
iteration. In one embodiment this information takes the form of a stored value for an 
additional bias. The additional bias is described in more detail below. In another 
embodiment this information takes the form of a stored value of the last correction 
term computed. The correction term is described in more detail below. 
[0099] In step 482, it is determined whether any of the other segments within 
the halo distance has had its bias changed in the last iteration. For example, it is 
determined whether any segment has a non-zero additional bias. If so, then the 
decision on whether the current segment is skipped or not depends upon whether the 
current segment was changed or skipped within either of the last two iterations. 
[0100] If the current segment has not been changed or skipped in the last two 
iterations, then the skip counter will equal 2. In step 484, it is determined whether 
the skip counter is equal to 2. If so, then it is assumed unlikely that the smaller 
changes in the last iteration on the neighboring segments will have a noticeable 
effect on the current segment, and the current segment is skipped. Control passes to 
step 485, where the skip counter is set back to 1 so that this current segment is not 
skipped again in the next iteration. Therefore, if the segment is skipped in this 
iteration, it will not be skipped in the next iteration. In the next iteration the current 
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segment will be visited again to determine whether its bias should be changed, as 
described in more detail below. If the current segment has not gone through two 
iterations without a change or a skip, the skip counter will be less than 2, and the 
segment will be processed in step 410 described below. 

[0101] If none of the other segments within the halo distance has had its bias 
changed in the last iteration, as determined in step 482 above, then control passes to 
step 4.86. In step 486, it is determined whether the bias of the current segment was 
changed in the last iteration. If so, it is assumed that the bias still needs some fine- 
tuning and the segment will not be skipped. In this case, the segment will be 
processed in step 410 described below. If it is determined in step 486 that the 
current segment was not changed during the last iteration, then there is no reason to 
expect the current segment will need further tuning. In this case, the current 
segment is skipped and control passes to step 438 described in more detail below. 

COMPUTING OPTIMAL BIAS WITH THE PROXIMITY EFFFECTS MODEL 
[0102 ] FIGs. 4E, 4F and 4G are flow diagrams showing details for 
embodiment 410a of the optimal bias determination step 410 in FIG. 4A for the 
current segment during one iteration. 
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GENERATING THE INITIAL PROFILE OF MODEL AMPLITUDES 
[0103] In step 411, the proximity effects model is rim for the current 
segment with the current value of the optimal bias. On the first iteration (1=0), the 
current value of the optimal bias is given by the pre-bias, if any, set in step 405 in 
FIG. 4B. Otherwise the optimal bias is zero on the first iteration. On subsequent 
iterations, the current value of the optimal bias is that computed on the previous 
iteration. In the following the current optimal bias is represented by Bj-i, to indicate 
the bias was computed on the previous iteration for the second iteration (1=1) and . 
later. On the first iteration, where 1=0, Bo-i means the pre-bias, if any, or zero. 
Details of embodiment 411a of step 41 1 are shown in FIG. 5B. In embodiment 
411a, the proximity effects model used provides a pre-image value as well as an 
amplitude and threshold, 

[0104] In step 530, stored values for the current segment, comprising model 
pre-image values at the evaluation point or at any supplemental evaluation points or 
both, are retrieved. These stored values do not represent a segment displaced from 
the corresponding edge in the design layer by the current optimal bias, but, instead, 
represent a segment displaced by the optimal bias determined on the iteration before 
the previous iteration. On the first iteration, (1=0), there are no stored values of pre- 
image. 

[0105] The supplemental evaluation points are points associated with the 
segment that are displaced from the evaluation point, usually perpendicularly to the 
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segment. The supplemental evaluation points are used in some embodiments to 
compute the size and sign of experimental biases, as is described in more detail 
below. In the preferred embodiment, two supplemental evaluation points are used. 
[0106] In step 532, changes to the bias associated with the current segment 
and with other segments within the halo distance are retrieved. On the first iteration, 
this stored change in bias is just the pre-bias, if any, set in step 405 in FIG. 4B. If no 
pre-bias was set in step 405 for this segment, then the change in bias is zero on the 
first iteration. On the second iteration (1=1), the change in bias is the optimal bias 
set on the first iteration (1=0). On subsequent iterations (I>1), this stored change in 
bias is the difference between the optimal bias computed during the previous 
iteration and the value of the optimal bias computed on the iteration preceding the 
previous iteration. In this embodiment, the halo distance used to select neighboring 
segments is the halo distance of the current iteration computed according to 
Equation 2 in step 450. During the first iteration, the halo distance is h max . 

[0107] In step 534, the kernel functions are convolved with the changes in 
bias determined in step 532 for the selected segments. The result is a computed 
change in pre-image at the evaluation point or supplemental points or both. Then 
the new pre-image value is obtained by adding the change to the stored pre-image 
value or values, if any, in step 536. A new value or values of intensity and, hence, 
new values of amplitude and printed edge position are then based on the new pre- 
image in step 538. In this embodiment, the number of kernel functions used is the 
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number N] computed according to Equation 3 during step 450. During the first 
iteration, the number of kernel functions used is M. 

[0108] The product of step 538 is the initial predicted edge of step 41 1 in 
FIG. 4E along with an initial profile of amplitudes at the evaluation point or at the 
supplemental evaluation points, if any, or both. By definition, the predicted edge 
resulting from running the proximity effects model at the supplemental evaluation 
points, if any, is substantially the same as the predicted edge resulting from running 
the proximity effects model at the evaluation point. 

USING THE INITIAL PROFILE OF MODEL AMPLITUDES 
[0109] In step 4 1 2, the initial predicted edge and the initial profile of model 
amplitudes are used to compute a first correction term CD(0) that is related to the 
distance between the predicted printed edge and the target edge. In one 
embodiment, if a calibrated model is used that provides a spatial deviation, that 
deviation is used as the first correction term CD(0). In another embodiment, the 
initial correction term, CD(0), is set equal to the conventional bias. In another 
embodiment, CD(0) is calculated using equations presented below, as shown in 
steps 514 through 5 1 8 in FIG. 5 A. In other embodiments, other correction terms are 
used. 

[0110] According to embodiment 41 2a, an estimated correction term CD is 
computed based on the difference between the amplitude at the evaluation point and 
the amplitude threshold value. Recall that the amplitude can be the spatial deviation, 
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in which case the amplitude threshold is zero. Alternatively, the amplitude can be 
the intensity and the amplitude threshold is the constant or variable intensity 
threshold. Alternatively, the amplitude can be the difference between the intensity 
and the intensity threshold, in which case the amplitude threshold is a constant, zero. 
Any of these types of models can be used by embodiment 412a. For simplicity, the 
embodiment is described as if the amplitude, I, is optical intensity and the amplitude 
threshold, T, is the variable intensity threshold. 

[0111 ] In embodiment 4 1 2a, the calculation of CD includes the slope of the 
amplitude profile at the evaluation point according to Equation 4, below. The 
correction term, CD, is computed as follows: 

CD = (I e -T)/s ' (4) 
where I e is the amplitude at the evaluation point, T is the amplitude threshold value 

for that amplitude, and s is the slope of the amplitude profile at the evaluation point. 
Equation 4 is an estimate of the distance from the evaluation point to the point where 
the amplitude profile attains the value of the amplitude threshold. Equation 4 is a 
very good estimate when the evaluation point is in a portion of the profile that is 
near linear, or the threshold crossing point is nearby. Equation 4 is a poor estimate 
when the evaluation point is in a curved portion of the profile, or far from the point 
where the profile attains the value of the threshold. For example, in FIG. 3B, 
Equation 4 would be a poor estimate of the correction needed on curve 330 if the 
evaluation point were at 490 nm; but Equation 4 provides a substantially better 
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of the profile near 400 nm. 

[0112] Different embodiments utilize different steps to obtain a value of the 
amplitude at the evaluation point Ig. In one embodiment, Ie is obtained directly by 
running the proximity effects model at the evaluation point. 

[0113] In an alternative embodiment, two supplemental evaluation points are 
defined, and the proximity effects model is run at the two supplemental evaluation 
points. The first supplemental evaluation point is positioned inside the segment a 
distance gi nn er from the evaluation point and has an amplitude of Ij . The second 
supplemental evaluation point is positioned outside the segment a distance g ou ter 
from the evaluation point and has an amplitude of I2. For example, in FIG. 3 A, the 
first supplemental evaluation point is placed at point 309d and the second 
supplemental evaluation point is placed at point 309c. In this case, I e is computed 
using.Equation 5 a. 

Ie = ( I| * gouter + h * ginner) / (ginner + gouter) (5a) 

[0114 ] In one embodiment, ginner is substantially equal to a specification 
tolerance for under-printing and g oute r is substantially equal to a specification 
tolerance for over-printing. An advantage of this embodiment is that a clear 
indication is provided when the predicted edge is within tolerance of the original 
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fabrication layout edge. That is, if the value of CD calculated according to Equation 
4 lies in the range from -gj nn er to gouter, then the predicted edge is within the 
specification tolerances for both over-printing and under-printing, without 
approximation. 

[0115] Where the proximity model utilizes look-up tables to implement 
convolutions of the kernel functions, the entries in the look up tables often 
correspond to spatial locations of edges separated by a model grid distance. 
Therefore, in one embodiment, gi nn er and g ou ter are selected so the two supplemental 
evaluation points are spaced apart by the model grid distance. In some embodiments 
ginner and go Ut er are selected to be substantially equal. 

[0116] Equation 5a can also be used with embodiments using the evaluation 
point with only one supplemental evaluation point outside the edge by setting ginner 
to zero and setting g ou ter to a non-zero positive distance and using the evaluation 
point amplitude for Ij . Similarly, Equation 5a can also be used with embodiments 
defining only one supplemental evaluation point inside the edge by setting gouter to 
zero and setting gj nner to a non-zero positive distance and using the evaluation point 
amplitude for I2. In either case, the non-zero positive distance is set substantially 

equal to the model grid distance in some embodiments. 
[0117] To obtain a value of the slope, s, at the evaluation point in the 
embodiments using supplemental evaluation points, the model is run for each profile 
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at the evaluation point and a supplemental evaluation point, or at two supplemental 
evaluation points. In either case, the slope is computed as follows: 

S-(I 2 -Ii)/(ginner + gouter) (5b) 

Equation 5b is used with the evaluation point and a supplemental evaluation point 
outside the edge by setting gj nner to zero and using the evaluation point amplitude for 
ly. Equation 5b is used with the evaluation point and a supplemental evaluation point 
inside the edge by setting gouter to zero and using the evaluation point amplitude for 
I2. The choice of inside or outside supplemental evaluation point depends on the 

attention to be paid to over-printing versus under-printing. 
[0118] When the slope is computed according to Equation 5b, the sign of 
CD given by Equation 4 is a reliable indicator of the direction of a correction needed 
in a bias. For example, if this CD is positive, it means that the printed feature is 
printing too small, and that a positive correction to the bias is needed to compensate. 
For another example, in FIG. 3B, the supplemental evaluation point located outside 
the edge at a distance of about 390 nm, has a larger value of amplitude than the 
evaluation point. Consequently, the slope s. will be positive, and the value of CD on 
curve 330 will be negative, because the amplitude at the evaluation point Ie is less 
than the threshold value T being subtracted from it. This indicates the edge is 
printing too big. A negative change in the bias is indicated. A negative change in 
bias moves the edge in the fabrication layout inside the edge in the design layer. In 
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FIG. 3B, this direction is toward the right. Consequently, the negative value of CD 
indicates the segment should be moved to the right of its position for curve 330. 
[0119] FIG. 5 A shows more details of the techniques used in embodiments 
411b and 412a of steps 411 and 412, respectively, to utilize Equations 1 and 2. In 
step 510, on the first iteration, two supplemental evaluation points for the current 
segment are defined, one inward from the evaluation point for the segment by a 
inner distance gj nne r> and one outward from the evaluation point by an outer distance 
g outer . In step 512, the proximity effects model is run at the two supplemental 
evaluation points for the current value of bias to output a predicted edge position and 
a profile with two amplitudes, one amplitude at each of the supplemental evaluation 
points. In other embodiments, the evaluation point is used in lieu of one of the 
supplemental evaluation points. 

[0120] In step 514, the amplitudes at the supplemental evaluation points are 
used to compute the amplitude, Ie, at the evaluation point according to Equation 5a, 

and to compute the slope, s, according to Equation 5b: In other embodiments, the 
evaluation point is used in lieu of one of the supplemental evaluation points. In step 
516, the difference is computed between the amplitude threshold value, T, and the 
amplitude at the evaluation point, I e . In step 51 8, the correction term CD is 
computed according to Equation 4. In step 520, the pre-image values at the 
supplemental evaluation points are stored. Also stored, in this embodiment, is the 
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value of the correction term CD. 

[0121] , The use of a pre-image saves computation time because the new 
contributions to the amplitude output by the model at the evaluation or supplemental 
evaluation points are computed only for the changes since the last iteration, and the 
resulting changes in pre-image values are simply added to the pre-image values 
already stored. During the first iteration, the pre-image is empty and the 
computation of amplitudes at evaluation points must include all geometries within 
the original layout. The resulting pre-images are stored. In the second iteration, 
only the contributions to pre-image from the biases are computed and added to the 
stored pre-image values. In the next and subsequent iterations, only the differences 
in biases, but not the biases themselves, are used to compute changes in the pre- 
image values near the evaluation points. As the number of segments that have their 
biases changed tends to decrease with later iterations, the new contributions to the 
pre-image also tend to decrease. Therefore the computational cost of adjusting the 
layout dramatically decreases. 

[0122] In step 413, the first correction term CD(0) is compared to the 
specification tolerance. The correction term is considered an indication of the order 
of magnitude of the change needed in the current bias for the current segment. Thus, 
if the first correction term CD(0) is within the specification tolerance, then the bias 
for the current segment is converging on a value sufficiently precise to satisfy the 
specification requirements. If the first correction term CD(0) is not within the 
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specification tolerance, then a significant change in the bias for the current segment 
is indicated. In this case, control passes to step 414 in which the skip counter is reset 
to 0. The skip counter is used to count the number of iterations that the current 
segment has not had significant changes in its bias and has not been skipped already. 
When control passes to step 414, a significant change in the bias of the current 
segment is expected, and therefore, the skip counter is reset. 
[0123] If the first correction term CD(0) is within specification tolerance, 
then the current bias of the current segment is sufficiently precise to satisfy the 
specification, as stated above. When the values of the correction terms CD for all 
segments are within the specification tolerance, iterations are stopped. This is the 
termination criterion tested in step 440 of FIG. 4 A as described above. However, in 
this embodiment, the bias for the current segment will still be adjusted unless the 
first correction term CD(0) is also within a target tolerance smaller than the 
specification tolerance. Requiring the current segment to have a correction term 
within the smaller target tolerance provides some margin to compensate for the 
proximity effects that couple changes in bias for each segment with the changes in 
bias of neighboring segments. Thus, if the first correction term CD(0) is within 
specification tolerance, then control passes to step 415 wherein it is determined 
whether the first correction term CD(0) is within the target tolerance. If CD(0) is 
also within target tolerance, then no further bias changes are necessary and control 
flows to step 430. Otherwise, control passes to step 416. 
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[0124] In step 4 1 6, the slope of the initial profile, determined using Equation 
5b during the computation of the first correction term CD(0) in step 412, is tested. 
For the bias of the current segment to be considered correctable on this iteration, the 
slope must have large enough magnitude to attain the amplitude threshold value 
within a reasonable distance. The geometry of other features in the design layer 
limit how far an edge can be moved. If the magnitude of the slope is not large 
enough, the amplitude threshold value cannot be attained in the space provided. The 
slope at the nominal edge is a measure of the image contrast, especially when the 
amplitude is optical intensity. 

[0125] In step 41 6, the slope is compared to a minimum contrast slope, s m i n . 
If the magnitude of the slope exceeds or equals the minimum contrast slope, this 
segment is deemed correctable during this iteration and control passes to step 472. 
[0126] If the magnitude of the slope is less than the minimum contrast slope, 
then the bias for the segment is not corrected during this iteration. If this happens in 
the first iteration, then the segment is considered uncorrectable during any iteration, 
and the segment is dropped to avoid any further visits. This reasoning is 
implemented in step 417, which determines whether the current iteration is the first 
iteration. If not, the segment is simply skipped and control passes to step 438. If it 
is the first iteration, control passes to step 418 where the segment is flagged as 
uncorrectable, and then control passes to step 438. A segment flagged as 
uncorrectable is not selected as the current segment during step 406 of FIG. 4A in 
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later iterations. 

[0127] If the edge is deemed correctable in step 416 because the slope is 
greater than the minimum slope, s min , then control passes to step 472. In step 472 it 

is determined whether the current optimal bias equals or exceeds the corresponding 
maximum set for positive and negative biases during step 407 of FIG. 4B. If so, 
then control passes to step 474 to determine whether a correction on the order of CD 
would result in still exceeding the corresponding maximum. If so, then the segment 
is not corrected further and control passes to step 438 to check for the next segment. 
[0128] If the current optimal bias does not equal or exceed the corresponding 
maximum, or if the correction tends to reduce the bias below the maximum, then 
control passes to step 419 in FIG. 4F to determine a second experimental bias at 
which to run the model. 
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GENERATING THE SECOND PROFILE OF MODEL AMPLITUDES 
[0129] In step 41 9, a second bias is derived by adding to the current value of 
the optimal bias a first additional bias, d(0), that is a function of the first correction 
term, CD(0). 

[0130] In one embodiment, the additional bias has the same sign as the first 
correction term CD(0), and has a magnitude of the experimental additional bias. In 
this embodiment, the experimental additional bias is a function of the iteration 
number, I, computed according to Equation 1 during step 450 in FIG. 4A. During 
the first iteration, with iteration number I equal to zero, the experimental additional 
bias is the maximum experimental additional bias. Therefore, in this embodiment, 
the first additional bias is given by the following equation: 

d(0) = sign(CD(0))* |d,| (6a) 

where the pair of vertical bars indicate the absolute value of the quantity between the 
bars. Recall that I is the iteration number. 

[0131] In another embodiment, the additional bias is a function, fd, of the 
slope, s, and the correction term, CD. In this embodiment, a high contrast edge is 
defined as one having a slope greater than a predetermined value Shigh* Then, if the 
slope, s, is greater than Shigh, the additional bias is made proportional to the 
correction term, CD. The estimate of CD according to Equation 4 is a good 
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indication of the expected change in bias when the edge has high contrast. On the 
other hand, when the contrast is lower than s n i g h> the value of CD is more sensitive 

to small changes in amplitude and may yield an additional bias that is too large, 
especially if the model has some noise near the target edge. An additional bias that 
is too small is preferable to an additional bias that is too large. For example, an 
additional bias that is too large may cause changes in neighboring edges that 
contaminate their biases for many iterations. So, when the slope , s, is greater than 
Smin but less than Shi gn , the additional bias is made proportional to s rather than to 
CD. These relationships are captured by the following equation defining function 
fd: 

( kr*CD if s > s hi gh c 

fd(s,CD)= | k2 * s ifs high >s>s min . (6b) 

0 ifs min >s 

where ki and k 2 are preset constant coefficients. In an embodiment using Equation 
6b, the constants k\ and k 2 are set during step 400 of FIG. 4B. In this embodiment, 
then, the first additional bias is given by the following equation: 

d(0) = fd(s, CD(0)) (6c) 

[0132] In a third embodiment, the additional bias is set equal to the 
correction term CD(0) unless this additional bias is too large. To ensure this choice 
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does not yield an additional bias that is too large, the additional bias is selected as 
the minimum magnitude of the additional biases computed using CD(0), Equation 
6a and Equation 6c. That is, in this embodiment, the first additional bias is given by 
the following equation: 

d(0) = sign(CD(0)) * MIN [ |d,| , | f(s,CD(0))| ] (6d) 

where MIN[] indicates an operation that selects the quantity having the smallest 
value of two or more quantities listed within the brackets and separated by commas. 
[0133] The second experimental bias, D2, is then the sum of the optimal bias 
from the last iteration, B\.\, stored as the "current" optimal bias, and the first 
additional bias d(0) as given by the following equation: 

D2 = B M +d(0) (6e) 

[0134] In step 420, the proximity effects model is run for the segment 
displaced from the corresponding segment in the original fabrication layout by the 
second experimental bias, D2. The proximity effects model is run at the evaluation 
point or the supplemental evaluation points, if any, or both, to output a second 
printed edge position and two amplitudes that constitutes a second profile of model 
amplitudes. Details of running the proximity effects model are similar to those 
described in FIG. 5B for step 411. Control than passes to step 42 1 . 
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USING THE SECOND PROFILE OF MODEL AMPLITUDES 
[0135] In step 421, the second profile of model amplitudes is used to 
compute a second correction term CD(1) using Equation 4 and Equation 5 a and 
Equation 5b (as shown in steps 514 through 518 in FIG. 5 A). 
[0136] In step 422, the signs of the first and second correction terms, CD(0) 
and CD(1), are compared. If the signs are different, then control passes to step 423 
where an interpolation is performed, as described below. If the signs are the same, 
then extrapolation is indicated. In this case, control passes to step 424 where 
extrapolation is avoided by using the first additional bias d(0) computed in step 419 
as the correction to the current optimal bias. Essentially, the first additional bias 
d(0) computed in step 419 is a maximum change in bias for this iteration. The new 
optimal bias for this iteration, Bi, is obtained by adding the first additional bias d(0) 
to the current optimal bias as was done in equation 6e, above. So, when 
extrapolation is indicated, the following equation used: 

B, = D2 = B M +d(0) (7a) 

[0137] In step 423, a linear interpolation is performed to obtain an optimal 
bias associated with a correction term equal to zero, where the amplitude at the 
evaluation point equals the threshold value. In this embodiment, the interpolated 
bias is a third bias computed as follows: 

D3 = Bj.| + d(0) * CD(0) / (CD(0) - CD(1)) (7b) 
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where CD(O) is the first correction term, CD(1) is the second correction term, d(0) is 
the first additional bias, and D3 is the third bias. The new optimal bias is based on 
D3. D3 is used as the new optimal bias if CD(0) and CD(1) have opposite signs, as 
shown in the following equation: 

B! = D3 (7c) 

[0138] In another embodiment, Equation 7b is also used in step 424 even if 
the two correction terms have the same sign and a maximum bias is set. Then, if the 
magnitude of D3 is less than the maximum bias, D max , set in step 407 for a bias of 
the corresponding sign, the optimal bias is set to D3. Otherwise, the optimal bias is 
set to the maximum bias, D max . This embodiment is given by the following 
equation: 

B, = sign(CD(0)) * MIN [|D3j, |D max | ] (7d) 
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THE TWO-PASS OPTION 

[0139] After the new optimal bias, Bi, is determined in step 423 or 424, it 
can be used in step 430 to displace the current segment using a one-pass option. 
Using a two-pass option, a third profile and second interpolation is attempted before 
control is passed to step 430. Therefore, in the illustrated embodiment, control 
passes to step 425 to determine whether the two-pass option has been selected. If the 
two-pass option has not been selected, control passes to step 430. 
[0140] If the two-pass option has been selected, control passes to step 426 to 
determine whether the second correction term CD(1) lies within the target tolerance. 
If so, there is little to be gained by computing the third profile and the second 
interpolation; so, control again passes to step 430. If the second correction term 
CD(1) does not lie within the target tolerance, then control passes to step 427 in FIG. 
4G to continue with the two-pass option. 

GENERATING THE THIRD PROFILE OF MODEL AMPLITUDES 
[0141] In step 427 of FIG. 4G, a further experimental bias is calculated 
based on the newly computed optimal bias, B\. In step 428, the proximity effects 

model is run for the segment displaced from the corresponding segment in the 
original fabrication layout by the further experimental bias. The proximity effects 
model is run at the evaluation point or the supplemental evaluation points or both to 
output a third predicted edge and the two amplitudes that constitutes a third profile 
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of model amplitudes. Details of running the proximity effects model are similar to 
those described in FIG. 5B for step 411. Control than passes to step 429. 

USING THE THIRD PROFILE OF MODEL AMPLITUDES 
[0142] In step 429, the third profile of model amplitudes is used to compute 
a third correction term CD(2) using Equation 4 and Equation 5 (as shown in steps 
514 through 518 in FIG. 5A). 

[0143] In step 461 , the sign of the third correction term CD(2) is compared 
to the signs of the first two correction terms, CD(0) and CD(1). If the signs are 
different, for example if either CD(0) or CD(1) has a different sign than CD(2), then 
a second interpolation is worthwhile. In this case, control passes to step 462 where a 
linear interpolation is performed between CD(2) and the closest correction term of 
opposite sign, to obtain the bias associated with a correction term of zero. If the 
signs are the same, then extrapolation is indicated. Control passes to step 463 where 
extrapolation is avoided by using the new optimal bias as computed in step 424. In 
this case, extrapolation was indicated in step 422 so that step 424 must have been 
performed. In either case, control passes to step 430 to displace the current segment 
by the new bias. 

DISPLACING THE CURRENT SEGMENT 

[0144] FIG. 4H is a flow diagram showing details for embodiment 430a of 
the displace segment step 430 in FIG. 4A. Step 43 1 denotes the relationship 
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between the current optimal bias B\. u the additional bias di, and the new optimal 
bias, Bj. That relationship is indicated by the following equation: 

di = B,-B,., (8) 
[0145] Step 43 1 need not involve a subtraction if the first additional bias 
d(0) computed in step 419 or step 427 or both is used and is still available, such as in 
memory or on nonvolatile storage. In that case, d\ = d(0). In step 432, a computer 

variable B c holding the current optimal bias for the current segment is set equal to 
the new optimal bias, Bi. In step 434, it is determined whether the additional bias, 
di, marked in step 431 is equal to zero. If so, control passes to step 436 where a step 
counter is incremented. Incrementing the step counter indicates that the current 
segment did not have its bias changed during this iteration. That is, the additional 
bias di, denoted in step 43 1 was zero. When the skip counter reaches 2, the segment 

is skipped in the next iteration, during step 408 of FIG. 4A. In any case, control 
passes to step 435. 

[ 0 14 6 ]. In step 435, it is determined whether the new optimal bias exceeds 
the corresponding maximum bias, D max , for the same sign. If so, then, in step 326, 

the magnitude of the new optimal bias is set to the corresponding maximum, D max , 

redefining B|, and the additional bias, dj is recomputed using Equation 8. 

[0147] In either case, control next passes to step 438 where the next segment 
in the layout is determined. 
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REPLACING STORED VALUES 

[0148] Only the latest values of optimal bias, pre-image, additional bias, and 
last computed correction term, CD, are stored at the end of each iteration. When the 
stored values are replaced, the steps that check neighboring segments for current 
bias and changes are presented with the replaced values. Thus time of replacement 
determines when corrections are felt by neighboring segments. 
[0149] In some embodiments, the stored values of optimal bias, pre-image, 
additional bias, and last computed correction term, CD, are replaced during step 430 
before going to the next segment. In other embodiments, the stored values are 
replaced after step 438 determines that there are no more segments in the layout to 
be visited on the current iteration. The timing of this replacement determines when 
corrections made on one segment are detected by other segments. 

SEPARATE TREATMENT OF CORNER SEGMENTS 

[0150] In another embodiment, corner segments including only one vertex of 
an edge in the original fabrication layout, are treated differently than non-corner 
segments. The optimal bias computed for a corner segment is significantly reduced 
if all the non-corner segments on a polygon of the original fabrication layout are 
corrected with their optimal biases before computing and applying the bias for a 
corner segment. 

[0151] For example, in FIG. 9, a polygon 910 in an original fabrication 
layout has four edges defined by vertices 912, 914, 916 and 918. The polygon has 
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been dissected into eight segments with dissection points that include the four 
vertices as well as non-corner dissection points 920, 922, 924 and 926. Each of the 
resulting eight segments 930, 932, 934, 936, 938, 940, 942 and 944 has an 
evaluation point, 931, 933, 935, 937, 941, 943 and 945, respectively. The proximity 
effects model predicts the polygon 910 will print the elliptical feature 950. FIG. 9 
depicts an optimal bias 960 computed for segment 938 at evaluation point 939 that 
moves that segment to position 964 in an adjusted fabrication layout after one 
iteration. According to some embodiments, the bias computed for corner segment 
940 including vertex 918 with evaluation point 941 will be based on a correction 
term CD proportional to the distance marked by arrow 963. 
[0152] According to another embodiment, a corner segment having one 
vertex, like corner segment 940, does not have its bias determined until after all the 
non-corner segments have had their biases determined and applied. Segment 938 
has two vertices and therefore is a line end segment, not a corner segment, according * 
to this definition of a corner segment. Therefore the bias 960 computed for segment 
938 is applied before computing a bias for corner segment 940. Segments 930, 934 
and 942 are also non-corner segments. To illustrate the effect of treating corner 
segments separately, it is assumed that a negative bias (not shown) is computed for 
non-corner segments 934 and 942 as well as the positive bias 960 computed for 
segment 938. After these biases are applied to the corresponding segments, the 
proximity effects model predicts the different elliptical feature 970 will be printed. 



60 



NTI-019-5-1D PATENT 



If the corner segment 940 is now treated to determine a bias at evaluation point 941, 
it will be based on a correction term having a size illustrated by arrow 973. As can 
be seen, the correction term 973 for the corner segment computed in this way is 
substantially smaller than the correction term 963 computed in the alternative 
embodiment. Thus, the fabrication layout can be adjusted using less severe biases in 
corner segments. 

[0153] According to one embodiment for treating corner segments ' 
separately, in step 406 of FIG. 4A, the next segment selected will be only non-corner 
segments on one iteration and only corner segments on the next iteration. In this 
embodiment, each current bias (B c ), pre-image values at the evaluation point or 
supplemental evaluation points or both, correction term (CD), and additional bias 
(d), replaces the stored values of those quantities for the segment in step 440 in some 
embodiments, and in step 450 in other embodiments. 

[0154] According to another embodiment for treating corner segments 

separately, in step 406 of FIG. 4 A, the next segment is chosen so that all non-corner 

segments of a polygon are visited before any corner segment of the polygon. In this 

embodiment, each bias, pre-image, correction term and additional bias replaces the 

stored values of those quantities for the segment when there are no more non-corner 

segments on the polygon, and the next segment on the polygon is a corner segment. 

* 

This check and replacement is made in step 438 in some embodiments, in step 406 
in other embodiments. 
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CORRECTION DAMPING 

[0155] Applying a bias that is too large to one segment can adversely impact 
other segments. This impact is not considered during the iteration in which the 
correction is made because the neighboring segments are corrected based on the 
stored prior value of the bias which is usually not updated until the end of the 
iteration. On the next iteration the affected segments sometimes have to be 
overcorrected to offset the impact. This results in oscillating over-corrections that 
can seriously impair convergence of the adjusted fabrication layout to a stable 
solution. This undetected coupling is particularly strong in the first few iterations 
where the additional experimental bias, dj, is large (see Equation 1). To reduce the 
effects of undetected coupling between neighboring segments, the computed 
correction is damped by multiplying by a damping factor p less than or equal to one. 
[0156] In one embodiment, the damping factor increases from a value of po 
toward a value of 1.0 as the iteration number I increases. The change of damping 
factor p with iteration is given by the following equation: 

pi= \-{\-ptffA x (9a) 

where I is the iteration number; p\ is the damping factor on iteration I; po is the 
minimum damping factor; and f4 is a fourth fading factor having a value greater 
than 1 . A value for the fourth fading factor is chosen in order to control the speed 
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with which the damping factor increases to a value of 1 .0. The larger the fourth 
fading factor, the greater the speed with which the damping factor increases. In this 
embodiment, values for po and f4 are set in step 402 of FIG. 4 A. 

[0157] Whenever a bias is computed using interpolation, the calculated bias 
is reduced based on the damping factor according to the following equation: 

B,= A*D3 (9b) 

where Bj is the new optimal bias for the current iteration and D3 is the third bias 
obtained via interpolation as given in Equation 7b. According to this embodiment, 
Equation 9b is substituted for Equation 7c in steps 423 and 462. 
[0158] In another embodiment, correction damping is achieved by 
interpolating to a bias associated with a CP value other than zero. In Equation 7b 
above, a value for a third bias D3 is computed by interpolating the additional bias 
between CD(0) and CD(1) to where CD equals zero. If there is concern that the 
optimal bias computed based on this D3 is too large, then a smaller D3 can be 
obtained by interpolating to a target value of CD closer to CD(0) than zero is to 
CD(0). This new target CD is defined.by a second damping factor y. That is, the 
new target CD for the interpolation is equal to the quantity y* CD(0) rather than 0.0. 
[0159] In one embodiment, the second damping factor decreases from a 
value of yq, less than one, toward a value of 0.0 as the iteration number I increases. 
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The change of damping factor ywith iteration is given by the following equation: 

n= (roV® 1 (ioa) 

where I is the iteration number; y\ is the second damping factor on iteration I; , 

less than 1.0, is the maximum second^damping factor; and f5 is a fifth fading factor 
having a value greater than 1. A value for the fifth fading factor is chosen in order 
to control the speed with which the second damping factor decreases to a value of 
0.0. The larger the fifth fading factor, the greater the speed with which the damping 
factor decreases. In this embodiment, values for /q and f5 are set in step 402 of FIG. 
4A. 

[0160] According to this embodiment, the interpolation described by 
Equation 7b is replaced by the following interpolation equation: 

. D3 = B M + d(0) * (1- n )♦ CD(0) / (CD(0) - CD(1)) m (10b) 

where the terms other than ^ are as described above with respect to Equation 7b. As 
the iteration number I increases and % decreases to 0.0, Equation 10b approaches 
Equation 7b. 

CONSIDERING NEIGHBORING SEGMENT BIAS CHANGES 

[0161] In many of the above embodiments, estimates of bias changes in 

neighboring segments are not considered until after the stored values are updated at 
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the end of each iteration. As mentioned above, this leads to a failure to account for 
coupling and can cause overly large bias changes that hamper convergence to a 
stable adjusted fabrication layout. 

[0162] According to another set of embodiments, estimates of changes in 
neighboring segments are considered. In these embodiments, steps 41 1 and 412 of 
FIG. 4E and step 419 of FIG. 4F are performed not just for the current segment but 
also for one or more segments closest to the current segment. In one embodiment, 
the non-corner segments on either side of the current segment and from the same 
edge as the current segment are considered together with the current segment. This 
way, the second predicted edge produced during step 420 for the current segment 
includes some of the coupling effects resulting from estimated changes in bias 
applied to the neighboring segments in step 419. After step 419, the second 
correction term CD(1) computed in step 421 and the new optimal bias computed in 
step 423 need only be performed for the current segment in this embodiment. 
[0163] Using these embodiments, when one segment is corrected, the 
potential changes in the environment of the segment are being taken into account. 
Thus some of the coupling is expected to be captured in the new optimal bias. By 
capturing more of the coupling in each iteration, the total number of iterations to 
achieve the tolerance specification is expected to decrease. Thus these embodiments 
are expected to speed the convergence to an optimal adjusted fabrication layout. 

AN INVERSE PROXIMITY EFFECTS MODEL 
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[0164] In another aspect of the present invention, the conventional proximity 
effects model is not used. Instead, an inverse proximity effects model is built and 
utilized to determine a bias for a segment in a fabrication layout. 
[0165] FIG. 6 is a flow diagram for building an inverse proximity effects 
model, process 600. The process includes step 602, in which a first set of edge 
positions are obtained for test patterns in a fabrication layout. In step 604, a second 
set of edge positions are obtained. The edge positions in the second set belong to the 
printed features in the printed features layer produced using a mask formed from the 
fabrication layout of step 602. It is not necessary in step 604 that the fabrication 
process be performed to obtain the second set of edges. Both the first set of edge 
positions and the second set of edge positions may be available already as a result of 
building the conventional proximity effects model for a particular suite of 
fabrication equipment. 

[0166] In step 606, kernel functions are fit to the first set of edge positions 
and second set of edge positions such that amplitudes related to positions in the first 
set of edge positions are output based on input of edges from the second set of edge 
positions. That is, the kernel functions are designed to output amplitudes related to 
edge positions in a fabrication layout when edge positions of the printed features 
layer are input. 

[0167] In step 608, an inverse proximity effects model is built by 
constructing a lookup table for each kernel function. A lookup table contains a 
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value for amplitude contribution at a position in the table related to a position in the 
printed features layer. Typically, the positions and orientations in the printed 
features layer are associated with rows and columns in the lookup table. The value 
in the lookup table at a particular row and column represents a convolution of the 
kernel function with an edge at a particular position in the printed features layer, and 
having a particular orientation. 

EXAMPLE EMBODIMENT WITH INVERSE PROXIMITY EFFECTS MODEL 
[0168] FIG. 7 is a flow diagram showing an iterative adjustment of a 
fabrication layout using an inverse proximity effects model, according to one 
embodiment of adjustment process 263 from FIG. 2. 

[0169] In step 702, parameters used by the techniques are initialized. In step 
704, the current iteration begins. A segment of the fabrication layout is visited only 
once per iteration. In step 706, one of the segments not yet visited in this iteration is 
made the current segment. In step 708, it is determined whether the current segment 
should be skipped. If the current segment is skipped, control passes to step 738 to 
determine whether another segment remains to be visited in the layout. If the 
current segment is not skipped, control passes to step 710 in which the inverse 
proximity effects model is run for the current segment. The inverse proximity effects 
model is run using the edges in the design layer as input. Then the output 
amplitudes indicating edges for the fabrication layout give the optimal positions to 
produce a printed features layer that agrees with the design layer. Also in step 710, 
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an optimal bias is determined based on the locations for the fabrication layouts 
indicated by the outputs from the model run. In step 730, the layout is adjusted if 
necessary by moving the current segment to a position separated from the 
corresponding segment in the design layer by the optimal bias. 
[0170] After it is determined in step 738, that another segment remains to be 
visited within the layout, control passes back to step 706 to make the next segment 
the current segment. If it is determined in step 738, that no other segment remains to 
be visited, then control passes to step 740 to determine whether all segments in the 
layout require corrections less than the specification tolerance. If so, then no further 
changes to the biases are required, and the process ends in step 790. If at least one 
segment requires a correction greater than the specification tolerance, then control 
passes to step 750 to increment to the next iteration, and visit at least one of the 
segments again. 

[0171] In another embodiment^ the inverse proximity effects model is not 
used iteratively on the layout, but is used once for each segment defined by an 
evaluation point and two dissection points in the layout. 

HARDWARE OVERVIEW 

[0172] FIG. 8 is a block diagram that illustrates a computer system 800 upon 
which an embodiment of the invention is implemented. Computer system 800 
includes a bus 802 or other communication mechanism for communicating 
information, and a processor 804 coupled with bus 802 for processing information. 
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Computer system 800 also includes a main memory 806, such as a random access 
memory (RAM) or other dynamic storage device, coupled to bus 802 for storing 
information and instructions to be executed by processor 804. Main memory 806 also 
may be used for storing temporary variables or other intermediate information during 
execution of instructions to be executed by processor 804. Computer system 800 
further includes a read only memory (ROM) 808 or other static storage device coupled 
to bus 802 for storing static information and instructions for processor 804. A storage 
device 810, such as a magnetic disk or optical disk, is provided and coupled to bus 802 
for storing information and instructions. 

[0173] Computer system 800 may be coupled via bus 802 to a display 8 1 2, 
such as a cathode ray tube (CRT), for displaying information to a computer user. An 
input device 814, including alphanumeric and other keys, is coupled to bus. 802 for 
communicating information and command selections to processor 804. Another type 
of user input device is cursor control 816, such as a mouse, a trackball, or cursor 
direction keys for communicating direction information and command selections to 
processor 804 and for controlling cursor movement on display 812. This input device 
typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second 
axis (e.g., y), that allows the device to specify positions in a plane. 
[0174] The invention is related to the use of computer system 800 for 
producing design layouts and fabrication layouts. According to one embodiment of 
the invention, optimal biases for fabrication layouts are. provided by computer 
system 800 in response to processor 804 executing, e.g., as threads, one or more 
sequences of one or more instructions contained in main memory 806. For example, 
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the optimal bias process runs as a thread 852 on processor 804 based on optimal bias 
process instructions 851 stored in main memory 806. Such instructions may be read 
into main memory 806 from another computer-readable medium, such as storage 
device 810. Execution of the sequences of instructions contained in main memory 
806 causes processor 804 to perform the process steps described herein. In 
alternative embodiments, hard-wired circuitry may be used in place of or in 
combination with software instructions to implement the invention. Thus, 
embodiments of the invention are not limited to any specific combination of 
hardware circuitry and software. 

[0175] The term "computer-readable medium" as used herein refers to any 
medium that participates in providing instructions to processor 804 for execution. 
Such a medium may take many forms, including but not limited to, non-volatile 
media, volatile media, and transmission media. Non- volatile media includes, for 
example, optical or magnetic disks, such as storage device 810. Volatile media 
includes dynamic memory, such as main memory 806. Transmission media includes 
coaxial cables, copper wire and fiber optics, including the wires that comprise bus 
802. Transmission media can also take the form of acoustic or light waves, such as 
those generated during radio-wave and infra-red data communications, 
[0176] Common forms of computer-readable media include, for example, a 
floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, 
a CD-ROM, any other optical medium, punchcards, papertape, any other physical 
medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, 
any other memory chip or cartridge, a carrier wave as described hereinafter, or any 
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other medium from which a computer can read. 

[0177] Various forms of computer readable media may be involved in carrying 
one or more sequences of one or more instructions to processor 804 for execution. For 
example, the instructions may initially be carried on a magnetic disk of a remote 
computer. The remote computer can load the instructions into its dynamic memory 
and send the instructions over a telephone line using a modem. A modem local to 
computer system 800 can receive the data on the telephone line and use an infra-red 
transmitter to convert the data to an infra-red signal. An infra-red detector can receive 
the data carried in the infra-red signal and appropriate circuitry can place the data on 
bus 802. Bus 802 carries the data to main memory 806, from which processor 804 
retrieves and executes the instructions. The instructions received by main memory 
806 may optionally be stored on storage device 810 either before or after execution by 
processor 804. 

[0178] Computer system 800 also includes a communication interface 818 
coupled to bus 802. Communication interface 818 provides a two-way data 
communication coupling to a network link 820 that is connected to a local network 
822. For example, communication interface 818 may be an integrated services 
digital network (ISDN) card or a modem to provide a data communication 
connection to a corresponding type of telephone line. As another example, 
communication interface 818 may be a local area network (LAN) card to provide a 
data communication connection to a compatible LAN. Wireless links may also be 
implemented. In any such implementation, communication interface 818 sends and 
receives electrical, electromagnetic or optical signals that carry digital data streams 
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representing various types of information. 

[0179] Network link 820 typically provides data communication through one 
or more networks to other data devices. For example, network link 820 may provide 
a connection through local network 822 to a host computer 824. Local network 822 
uses electrical, electromagnetic or optical signals that carry digital data streams. The 
signals through the various networks and the signals on network link 820 and 
through communication interface 818, which carry the digital data to and from 
computer system 800, are exemplary forms of carrier waves transporting the 
information. Computer system 800 can send messages and receive data, including 
program code, through the network(s), network link 820 and communication 
interface 818. 

[0180] The received code may be executed by processor 804 as it is 
received, and/or stored in storage device 810, or other non- volatile storage for later 
execution. In this manner, computer system 800 may obtain application code in the 
form of a carrier wave. . 

[0181] In the foregoing specification, the invention has been described with 
reference to specific embodiments thereof. It will, however, be evident that various 
modifications and changes may be made thereto without departing from the broader 
spirit arid scope of the invention. The specification and drawings are, accordingly, 
to be regarded in an illustrative rather than a restrictive sense. 
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